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The exact solution of A = 2 supersymmetric QCD found by Seiberg and 
Witten includes a prediction for all multi-instanton contributions to the low- 
energy dynamics. In the case where the matter hypermultiplets are massless, 
the leading non-perturbative contribution at weak coupling comes from the 
two-instanton sector. We calculate this contribution from first principles 
using the exact two-instanton solution of the self-dual Yang-Mills equation 
found by Atiyah, Drinfeld, Hitchin and Manin. We find exact agreement 
with the predictions of Seiberg and Witten for all numbers of flavours. We 
also confirm their predictions for the case of massive hypermultiplets. 



In recent years, the idea of electric-magnetic duality has lead to a dra- 
matic improvement in our understanding of gauge theories with extended 
super symmetry. Due to the work of Seiberg and Witten |], 0, and subse- 
quent generalizations, exact results are available both for the spectrum of 
BPS-saturated states and for the low-energy effective action. A common 
feature of these results is that they have highly non-trivial consequences at 
weak coupling which can be compared with the results of semiclassical cal- 
culations. In the case of N = 2 supersymmetric Yang-Mills (SYM) theory 
with gauge group SU(2), the exact low-energy effective action obtained in 
[p]] contains predictions for all multi-instanton contributions to low-energy 
Green functions. These predictions have been checked by explicit calculation 
in the one-|J and two-instanton jp sectors. In this letter we extend these 
calculations to the N = 2 theory coupled to Np hypermultiplets of matter in 
the fundamental representation. Our results provide a semiclassical test of 
the exact solution of this model at low-energy proposed in 0. The analysis 
relies heavily on the results for SYM theory (ie Np = 0) derived in || and we 
will concentrate primarily on the new features which appear in the presence 
of matter hypermultiplets. A more complete account of these calculations 
will appear elsewhere [[|. 

We begin by restricting our attention to the case of massless hypermul- 
tiplets. In this discrete symmetry ensures that only even numbers 
of instantons contribute ||. For Np < 4, where the theory is asymptoti- 
cally free, the exact low-energy effective action, including all terms with at 
most two derivatives or four fermions, can be derived from a holomorphic 
prepotential which has the following expansion at weak coupling ||, 

where \1/ is an N = 2 chiral superfield containing the massless photon and its 
superpartners, Ajv f is the dynamically generated scale of the theoryQ and C 

1 The numerical values of the constants F^k depend on the definition of the scale 
Ajy F - In this paper, as in we adopt the A-parameter of the Pauli-Villars regularization 
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is a numerical constant. The logarithmic term comes from one-loop perturba- 
tion theory while the k-th term in the infinite series corresponds to the leading 
semiclassical contribution of 2k instantons. In the following we will only con- 
sider the first term in this series: the two-instanton contribution. Including 
the result for Np = 0,0 a straightforward Taylor expansion of the exact so- 
lution specified in 0, yields the numerical values J 7 ^ = 5/2 4 , J 7 ^ = — 3/2 5 
and = 1/2 6 . The coefficient T^f 1 corresponds to an additive constant 
in the prepotential which does not contribute to the low-energy effective La- 
grangian and is irrelevant for our purposes. Following |7|, [| |4[] , we will check 
these values by comparing the corresponding contributions to the four anti- 
fermion Green function G^ F (xi,X2, ^3,^4) = (A^i)^^^^)^^)) with 
the result of an explicit two-instanton calculation. 

We also consider the case Np = 4, where the perturbative /3-function 
vanishes. In this case the theory is believed to have an exact electric-magnetic 
duality which requires that all non-perturbative corrections to the low-energy 
effective action should also vanish p[. In the following, we explicitly check the 
vanishing of the two-instanton contribution to Gf\ Finally we will briefly 
describe how these results are modified by the introduction of masses for the 
hypermultiplets. 

In the absence of matter, the particle content of the N = 2 theory is con- 
veniently described in terms of two N = 1 superfields which transform in the 
adjoint representation of the gauge group. The gauge superfield V contains 
the Yang-Mills gauge field v m and the gaugino A, while the complex chiral 
multiplet $ contains the adjoint Higgs A and the Higgsino ip. The theory has 
flat directions along which the scalar field A acquires a vacuum expectation 
value (VEV), A d , breaking the gauge group down to an abelian subgroup: 
SU{2) — > U(l). Without loss of generality we may choose A c \ = vt 3 /2 where 
v is an arbitrary complex number which parametrizes the moduli space of 
gauge-inequivalent vacua. In addition we introduce Np matter hypermul- 
tiplets which transform in the fundamental representation of 577(2). Each 

scheme which is appropriate for instanton calculations. This is related to the A-parameter 
of !|a S 4A^ = (A^) 4 -^ |. 

2 Note that, for Np = 0, the prepotential has contributions from all numbers of 
instantons. 
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N = 2 hypermultiplet corresponds to a pair of N = 1 chiral multiplets, Qi 
and Qi where i = 1,2 . . . Nf, which contain scalar quarks (squarks) and g« 
respectively and fermionic partners Xi an d Xi- Importantly we are restricting 
our attention to the Coulomb branch of the theory where the squarks do not 
acquire a VEV. In the N = 1 language, the matter fields couple to the gauge 
multiplet via a superpotential, 

N F ^ 

W = J2 V2Q t ^Qi + m.QiQ, (2) 

i=l 

suppressing colour indices. We have also included masses m 8 for the hyper- 
multiplets. Until further notice we set rrii = 0. 

As in [|J , we begin by constructing the short-distance superinstanton. The 
short- distance behaviour of the gauge field is determined by solving the self- 
dual Yang-Mills equation, while the fermion fields obey the covariant Dirac 
equation in the self-dual background. We must also solve the Euler-Lagrange 
equations for the scalar fields, subject to vacuum boundary conditions. In 
contrast, the large-distance behaviour of the massless fields is obtained by 
solving the linearized equations appropriate to the low-energy abelian theory. 
These equations yield power-law behaviour with arbitrary coefficients which 
must be fixed by matching onto the short- distance solutions || . As usual the 
resulting field configurations are not solutions of the exact coupled equations 
of motion but should be thought of as approximate solutions corresponding 
to quasi-flat directions in configuration space ||. 

The general solution of the self-dual Yang-Mills equation was found by 
Atiyah, Drinfeld, Hitchin and Manin (ADHM) |10| and is conveniently given 
in terms matrices with quaternionic entries Q. Specifically, the solution for 
the gauge field v m with topological charge n is given by, 

Vrn = U X d m U X U X U X = I (3) 

3 Quatcrnions are represented as 2 x 2 matrices; e.g., B aa = B m a%" , m — 0,...,3 
where ct^ q are the a-matrices of Wess and Bagger O] . Wherever there is no ambiguity 
we suppress these indices. For a full list of our conventions see Sec. 6 of B. 
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where / is the unit quaternion and U is an x-dependent (n + l)-component 
column vector with quaternionic entries which obeys, 

A fc A^A = A Afc = a Xk + b Xk x (4) 

where a and b are (n + 1) x n matrices of quaternionic constants. These 
constant matrices are further restricted by a nonlinear constraint which can 
be written as 

A kX A xl = {r l ) kl I (5) 

This equation constrains a and b by demanding that each entry in the matrix 
product on the LHS be proportional to the unit quaternion for all x. In gen- 
eral, after implementing the constraint and eliminating all redundant degrees 
of freedom the resulting solution depends on 8n independent parameters. In 
an asymptotic limit, these parameters can be identified as the positions, scale 
sizes and global gauge orientations of n well-separated BPST instantons. 



Specializing to the case of topological charge n = 2, the canonical form for 
the ADHM solution is specified by the following choices for a and b; 



(6) 



In general there are an infinite number of ways of solving the constraint (|5|). 
A particularly convenient choice is to eliminate the off-diagonal element a\ 
as p[; 

Q i = 7TZ i 2 Q 3(w 2 Wi - Wiw 2 ) (7) 





4|a 3 | 

In the limit 1 0,3 1 — > 00, the corresponding field configuration resembles two 
well-separated instantons centered at — (xq + a^) and — (xq — (13), with scale- 
sizes \wi\ and \w2 \ respectively. 



At short distance, the adjoint fermions A and ip take their zero modes values 
in the self-dual background. The complete set of 4n linearly independent 
solutions of the adjoint Dirac equation in the general ADHM background 
was found by Corrigan et al ||13|| . The general solution for each species of 
adjoint fermion depends on an (n + 1) x n matrix with constant Grassmann 
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spinor entries. Following Sec. 7.1 of M, we denote these matrices M 1 and 
A/" 7 for the gaugino A and the Higgsino if), respectively. The matrices A4 7 
and AC satisfy the following linear constraints, 



^"7 



M 7 



-M lT a« 



LW, = Af^b^ 



(8) 
(9) 



where we have suppressed the ADHM indices. After imposing the constraints 
we find a total of 8n independent Grassmann parameters. These parameters 
are fermionic collective coordinates which must be integrated over. In order 
to find the correct measure for this integration, it is necessary to evaluate the 
8nx8n normalization matrix for the adjoint zero modes. This is accomplished 
in Appendix B of M. 



Specializing to the case n = 2, the matrices M y and AC are parametrized 



as, 



Mr 



( /il 7 /i2 7 ^ 

4£ 7 + A^ 37 M l7 
V M ly 4£ 7 - M 3y J 



AC 



4e;+A4 7 ^v l7 



^2 7 



'1 7 "5 7 

This form already satisfies the constraint (||). Just like its bosonic counter- 
part in (P), the remaining fermionic constraint (§) is solved by eliminating 
the off-diagonal elements, M\ y and A/"i 7 ; 

1 



Mi 



2|a 3 | 2 
1 

2R 



a 3 (2aiA^ 3 + w 2 ^\ - Wiji 2 ) 
-a 3 (2aiM 3 + w 2 v\ - W\v-i) 



(10) 



The fermionic zero modes parametrized by £ Q and £' a correspond to the four 
supersymmetry generators which act non-trivially on the self-dual gauge field. 

In addition to the adjoint fermions which belong to the gauge multiplet, 
we must also consider the zero modes of the fundamental fermions x% an d 
Xi, for i = 1, . . . , Np, in the ADHM background. For arbitrary instanton 
number, these modes are given as ]14|, [15 



U x a b\kfkifcu 
U x a bxkfkifcu 



(11) 
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with a a Weyl and (3 an 577(2) colour index. The normalization matrix of 
these modes is diagonal in ADHM indices [15]. 

d 4 x(x?f(Xaj)p = ifKuKi, (12) 

In the two instanton case we have a total of ANp linearly independent zero 
modes, corresponding to the Grassmann parameters Ku and Ku with I = 
1,2. For later convenience we define the following, frequently occurring, 
combinations of the collective coordinates introduced above; 



L = 


\wi 


2 + 


\w 2 \ 


H = 


\wi 


2 + 


\w 2 \ 



+ 4|ai| 2 + 4|a 3 | 



Y 
Z 



Ml ^2 

N F 

i=l 



v x n 2 + 2MsNx - 2J\f 3 Mi 



(13) 



Following [^j we next solve the equation of motion for the adjoint scalar 
fields in the presence of the fermionic sources described above Q 

(v 2 a)\ = iv^KAy'a - ^y*) ^ 14 ) 

a Nf 

{v 2 rf)\ = -^£(x7x* + x7xtt) (is) 

displaying colour and flavour but suppressing Weyl indices. Note that, in 
the path integral, the fermion bilinear contributions to A and to are 
completely independent. The above equations must be solved subject to 
the vacuum boundary conditions A(x) — > A d = vr 3 /2 as \x\ — > oo. In [|], 
we described the method for solving these equations in detail for the case 
Nf = 0. The only difference in the present case is the presence of a source 
term on the RHS of fllSD involving the fundamental matter fermions. This 
term can be included simply by an appropriate redefinition of the constant 
parameters appearing in the solution for iV^ = 0. Explicitly we find, 

%A = U X (A + B) XX ,U X , (16) 
itf = U X (A + C) XX ,U X , (17) 

4 These equations follow directly from the Lagrangian (2.3) of Q with couplings to 
matter implied by the superpotential (Q). 
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Here A(v) = — A(v) and 



A\y = i5xoSyoA ci + kbxkekibix' (18) 
Bxx> = ^(MfM T -MfM T )xx> + Kbx k e kl b lxl (19) 
Cxw = pbxktkibix' (20) 



where 



with cu = W2A c \Wi — WiA c \w 2 - 

As mentioned above, these field configurations are not exact solutions to 
the coupled equations of motion. It follows that the Euclidean action evalu- 
ated on the configurations found above will depend explicitly on the collective 
coordinates. This dependence has two important effects. First it provides 
an effective IR cut-off on the integral over instanton scale sizes. Second, the 
dependence on fermionic collective coordinates lifts a subset of the fermion 
zero modes. For a SUSY gauge theory with an adjoint scalar, it turns out 
that all fermion zero modes are lifted except those which are protected by 
supersymmetry In the present for Np = 0, we will find that the 

only unlifted modes are the four modes parametrized by £ a and £,' a identified 
above. 



In ||, we obtained an exact expression for the super- instanton action for 
the case N F = 0. In that case, it was possible, using integration by parts 
and the equations of motion, to express the total action as a surface term. 
The resulting action was then determined soley by the 1/x 3 term in the 
expansion of adjoint scalar field at large distance. The final answer (Eqn (8.7) 
of Q]) contained a bosonic part and a part bilinear in the adjoint Grassmann 
parameters M.^ and M 1 . In the present case, the same method can be used 
to calculate the bosonic and fermion bilinear terms in the action. However 
there is also a new feature for iVV > 0: a term quadrilinear in the Grassmann 
parameters A / i 7 , A/" 7 , /Q and JCi, which cannot be written as a surface term in 
any obvious way. Instead this term is uniquely determined by the requirement 
that the total action should be a supersymmetric invariant 0. Including this 
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term, the super-instanton action is given as Sn f = 16n 2 /g 2 + Sn f where 

S Np = 4ir 2 L\v\ 2 - 16n 2 (u + XZ) (k + «) + AV2n 2 igu k A clf i k (22) 

with A = ig /8y/2. The only difference from the Np = expression is the 
dependence on the fundamental parameters /Q and Ki introduced by the 
replacement 5 uj — > to + XZ. 



At large-distances, \x\ 3> l/g\v\, the massless components of the anti- 
fermion have a power-law decay with a coefficient which can be related to 
the collective coordinate action Q2"i2] ) M [] 

_ ac 

^ D (x) = iV2^^ a£t (x,x ) (23) 

where S a a is the free propagator for a massless Weyl fermion. Here we are 
treating v and v as independent variables and using the fact that the action 
( |22|) is linear in each of these variables. Note that these fields now have a 
contribution trilinear in the Grassmann collective coordinates. 



The final ingredient in our calculation is the collective coordinate measure 
[ |16|j . An explicit expression in the case iVV = was given in Sec. 8.1 of 
Apart from the modification of the two-instanton action described above, the 
only new feature here is the integration over the parameters K,u and K,u. The 
accompanying Jacobian follows immediately from the normalization relation 
(|T2"|). For Np flavours we have, 

2 10 A 2(4 ~ iVF) 

/ dfi2 NpS) = 7{oT§ \ c / d 4 xod A a^d 4 ''Wid A W2 X d 2 ^d 2 M.j,d 2 \i\d 2 \i2 



5 In contrast it is easy to show that there is no dependence on the parameters ICi and 
Ki in the action of a single superinstanton. It follows that the Grassmann integrals over 
these parameters in the one-instanton measure are not saturated and the corresponding 
contribution to G$ vanishes as expected from (|l|). 

6 In particular, see the more recent, slightly expanded, version of Q where all expres- 
sions are given in the general case of complex v. 
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N F 

x d?£'d 2 N 3 d?v 1 d l V2 x X\_dK,iidK,2idK,iidlC2i 

i=i 



x 



exp (-Sjv f ) 



|a 3 | 2 - |ai| 2 | 
H 



where S2 = 16 is a symmetry factor associated with a discrete redundancy 
in the chosen parametrization of the two instanton solution [|J pl| . 

Putting the pieces together, we find for the four-point function; 

G { S(x 1 ,X2,x^x i ) = [ d f ii N ^f(x 1 )Xf(x 2 )^(x 3 )^(x i ) (24) 



The only remaining problem is to evaluate the (28 + 4iVV)-dimensiona][] in- 
tegral over the collective coordinates. The £ Q and £' a dependence in the four 



anti-fermion fields (p3|) saturates the integrals over these parameters in the 
measure. The remaining Grassmann integrals must be saturated by expand- 
ing the exponential of the two-instanton action or by the trilinear terms in 
(|23|). The bosonic integrals in (|24D can be performed using the same changes 
of variables used in the Np = case (see Sec. 8.2 of j4|). A straightforward 
calculation yields, 



Cat /A \ 2 ( 4 ~ N f) p 

G < $ F (x 1 ,x 2 ,x 3 ,x A ) = J d 4 x e a(3 S a a(xi,x )S l3 ^(x2,x ) 

x e 7<5 S 7i ,(x 3 ,x )S 5( 5(x4,xo) (25) 



where c = -945/8vr 2 , c x = 315/32tt 2 , c 2 = -15/64vr 2 and c 3 = c 4 = 0. We 
can now compare this result with the contribution to G$ from the pointlike 
four-fermion vertex in the exact low-energy effective Lagrangian 0, §]• Q ; 

(4) 1 d^J- r 

G n f (x\, X2, x%, X4) = J d^Xo e a ^S a a(xi, Xo)S /3 ^(x2, Xq) 

x e l5 S^(x 3 ,x )S S g(x 4 ,xo) (26) 

Substituting the expansion ([!]) for T in (|26D and comparing the two-instanton 
term with ( ^5|) we extract the values !F 2 = 5/2 4 , J 7 ^ = — 3/2 5 and = 
1/2 6 in precise agreement with the exact solution of Seiberg and Witten. 

7 We are omitting the integral over Xq in this counting. 
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In the exact solution, the two-instanton contribution to the four-point 
function vanishes automatically for Np = 3 and 4. This follows from the 
action of the four derivatives with respect to v in (2~B) on the corresponding 
terms in the prepotential (|l|). In fact, the vanishing of the semiclassical result 
(|25|) for Np = 3 and 4 has a similar explanation. Because of the asymptotic 
forms (T 



_J) of the anti-fermion fields, the semiclassical expression ( 24]) for 
Gtf a l so contains four derivatives with respect to v. By interchanging the 
order of these derivatives and the collective coordinate integrations in ( 24]) it 
is easy to show that C3 = C4 = in (E5|). 



In the remainder of this paper, we will briefly indicate how these results 
can be generalized to include non-zero masses, m^ F > mM F -i > ■ ■ ■ > mi, 
for the hypermultiplets. We begin by describing the proposed exact results 
which we wish to test. As in the massless case, the instanton expansion of the 
exact prepotential can be deduced from the elliptic curves given in (see 
also fll7] ). An important restriction on the form of the prepotential comes 
from the decoupling limit m^ F — > 00. In this limit the exact prepotential 
obeys; 

J**'>(v; W.a*,) ^jr(N F -i ){v . {mihANF i) (27) 

where the A-parameters for different numbers of flavours are related as, 

ms^N?' - A ^- F i (28) 

with A being the A-parameter of N = 2 SYM theory. Note that the relation 
( |27|) holds independently at each order in the instanton expansion. 



When all the masses are non-zero, the prepotential gets contributions from 
all numbers of instantons. In particular, there is a one-instanton contribution; 



^)(v; W,A^ 



. A 4-7V F 



^ (o) n 



Np 



n=l 



TT V^ 



(29) 



This clearly obeys the decoupling relation (27) and, in particular, the coeffi- 
cient JF^ -* can be extracted from the instanton expansion of the prepotential 
of the Np = theory: J 7 ^ = 1/2 M, |. 
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The exact two-instanton contribution can be expanded in terms of the 
following set of SO(2iVi?)-invariant polynomials in the masses, 



M, 



(N F ) 



Mj 



(N F ) 



M. 



(N F ) 



1 

Np 

E 

8=1 

Np 

El 2 
m i m j 

i<j 



M 



(Np) 
Np 



Explicitly we have, 
^W(v; { mi },A NF ] 



n=2 



7T 



Np 

n 



a v 2(4-7V F ) AT F 
5=0 



M. 



,2i) 



(30) 



(31) 



In the chiral limit, rrii — > 0, we recover the coefficients of (Q), hence we have 
/o = Jz^- The decoupling relation (|2T ) implies that the numerical 



coefficients are not independent but obey; 



(JV F ) 



(Np-1) 



J 6-1 



Jo 



(Np-5) 



(32) 



It follows that the constant f$ Np ^ is equal to the coefficient J 7 ^* ^ °f t ne 
massless case. Hence, to extend our check of the two-instanton prediction to 
the massive case, it suffices to verify ([H]) and (|32"|). 



The inclusion of a mass term in the instanton calculation is straightfor- 
ward. For any number of instantons, the only modification is to introduce 
an additional dependence on the Grassmann parameters Ku and JCu in the 
instanton action. Specifically, the mass term for the fundamental fermions 



X and x can De evaluated on the zero modes ( |TT|) using fll2D and gives the 
additional term; 

N F 



7T 



(33) 



i=i 
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The expressions (p9|), fl3~T|) and fl32| ) can be verified by calculating the corre- 
sponding one- and two-instanton contributions to G$ as before. The only 
modification to the massless calculation comes from expanding the exponen- 
tial exp(— S mass ) prior to performing the integrals over JCn and Ku. More 
complete details of all the calculations described above will appear in the 
near future . 

ND is supported by a PPARC Advanced Research Fellowship, ND and 
VVK are supported by the Nuffield foundation. MPM is supported by the 
Department of Energy. 
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